Tagging the chemical evolution history of the Large Magellanic 

Cloud disk ^ 



. . . Emilio Lapenna 
(N" 

^1 Dipartimento di Astronomia, Universitd degli Studi di Bologna, Via Ranzani, 1 - 40127 

(N : Bologna, ITALY 
>[ 

0| eniilio.lapenna2@unibo.it 

CN ■ Alessio Mucciarelli 

Dipartimento di Astronomia, Universitd degli Studi di Bologna, Via Ranzani, 1 - 40127 

<d ■ Bologna, ITALY 

O' 

^' 
Of 

Q . Livia Origlia 

in: 

■ INAF - Osservatorio Astronomico di Bologna, Via Ranzani, 1 - 40127 Bologna, ITALY 

Cd ■ 

livia.orifflia@oabo.inaf.it 

>: 

O . and 

o 



alessio . mucciarelli2@unibo . it 



Francesco R. Ferraro 



Dipartimento di Astronomia, Universitd degli Studi di Bologna, Via Ranzani, 1 - 40127 
■ Bologna, ITALY 

f rancesco . f erraroSSunibo . it 
5^ . ABSTRACT 



We have used high-resolution spectra obtained with the multifiber facihty 
FLAMES at the Very Large Telescope of the European Southern Observatory to 
derive kinematic properties and chemical abundances of Fe, O, Mg and Si for 89 
stars in the disk of the Large Magellanic Cloud (LMC). The derived metallicity 
and [a/Fe], obtained as the average of O, Mg and Si abundances, allow us to 
draw a preliminary scheme of the star formation history occurred in this region 
of the LMC. The derived metallicity distribution shows two main components: 
one component (comprising ~ 84% of the sample) is peaked at [Fe/H] = -0.48 
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dex and it shows an [a/Fe] ratio slightly under solar ([a/Fe] ~ -0.1 dex). This 
population was probably originated by the main star formation event occurred 
3-4 Gyr ago (possibly triggered by tidal capture of the Small Magellanic Cloud). 
The other component (comprising ~ 16% of the sample) is peaked at [Fe/H] ~ -1 
dex and it shows an [a/Fe] ~ 0.2 dex. This population was probably generated 
during the long quiescent epoch of star formation in between the first episode and 
the most recent bursts. Indeed, in our sample we do not find stars with chemical 
properties similar to the old LMC globular clusters nor to the iron-rich and a- 
poor stars recently found in the LMC globular cluster NGC 1718 and predicted 
to be also in the LMC field, thus suggesting that both these components are 
small (< 1%) in the LMC disk population. 

Subject headings: stars: abundances — Magellanic Clouds — techniques: spec- 
troscopic 



Introduction 



Despite several years of studies, the chemical evolution history of the Large Magellanic 
Cloud (LMC) as well its star formation history (SFH) is still poorly understood. The LMC 
has experienced a complex SFH due t o the interaction occu rred both with the Small Mag- 
ellanic Cloud (SMC) and the Galaxy (IBekki fc Chiball2005l ). A clear evidence of this quite 
complex evolution can be recognized in the LMC star cluster system, characterized by a wide 
range of age and metallicities. The study of these clus ter stellar populatio i is reveals the exis- 
tence of at least three components: the old (^ 13 Gyr,| 

120061 . 



metal-poor component ([Fe/H] < -1.5 dex, [Olszewski et al. 



1991 



iGrocholski et al. 



Mucciarelli et al.l 120101 ) probably formed d uring the first episode of star forrnation (SF), the 



intermediate age component (~ 1-3 Gyr, iGallart et al.ll2003l. 



the d ominant one, and a young component (< 1 Gyr, 
20 111 ) that include the last formed clusters. 



Brocato et al 



2003 



Ferraro et al. 2004) which is 



Mucciarelli et al. 



The SFH of field stellar populations is less known. ISmecker-Hane et al.l (|2002| ) found 
that the dominant stellar population in the central bar was formed betwe en 4-6 Gyr and 
1-2 Gyr ago. This result was substantially confirmed by the simulations of iBekki &: Chiba 



^Based on observations collected at the ESO-VLT at Cerro Paranal (Chile) under the program 080. D- 
0368(A). 
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( 120051 ) who found that the formation of young stellar populations in the LMC bar is asso- 
ciated with efficient SF in the last few Gyr (in particular, ^ 2 Gyr a g o). A comprehensive 
spectroscopic study of field populations was performed by ICole et al.l (|2005[ ) measuring the 
infrared Call triplet of 373 giant stars located around the center of the bar. They found a 
distr ibution peaked at [Fe /H] ~ -0.4 dex with a low metallicity tail reaching [Fe/H] < -2.0 
dex. ICarrera et al.l (120081 ) have measured the infrared Call triplet in four fields at different 
radial distances (~ 3°, 5°, 6° and 8°, that is between 2.6 and 7 kpc) from the center of 
LMC. They found an average [Fe/H] ~ -0.5 dex for the closest fields, a lower [Fe/H] ~ -0.8 
dex was found only in the more distant one. A detailed spectroscopic anal ysis of 59 giant 
stars l ocated in the inner disk at ~ 1.2 kpc from the center was performed by lPompeia et al. 
( 120081 ) ■ who found a distribution peaked at [Fe/H] ~ -0.75 dex. 

This paper is part of a project devoted to investigate the kinematic and chemical proper- 
ties of the stellar populations in the LMC through the use of high-resolution spectra, able to 
provide accurate information about the kinematics, the metallicity and the chemical abun- 
dance of individual el ements of these stars. Previo us papers of the pr oject have discussed 



the properties of old (IMucciarelli et al. 



2009 



20101) . interrnediat e-age ( iFerraro et al.l 12006 



Mucciarelh et al. 2008 ) and voung (jMucciarelli et al. 2011, 2012) globular clusters (GC) in 



the LMC. This paper is the first of the project dedicated to the kinematic and chemical 
characterization of the LMC field stellar populations: we present chemical patterns for a 
sample of 89 Red Giant Branch (RGB) stars members of the LMC and located in the field 
around the old metal-poor GC NGC 1786. This represents to date the largest sample of field 
giants in the LMC for which high-resolution spectra have been obtained. 



Observations 



We have observed 91 stars located i n the region surr ounding the old GC NGC 1786 at 
~ 2° NW from the center of the LMC (iKim et al.lll998). The sp ectra have been acquired 
with the multi-object spectrograph FLAMES ( jPasquini et al.ll2002l ) at the Kueyen ESO-VLT 
telescope. 

The spectroscopic targets have been selected by using the near-infrared (J, H and K 
bands) photometric catalog, obtained by combining the SOFI catal og for the inner ~ 2.5 
arcmin, corresponding to the area covered by the GC NGC 1786 (see IMucciarelli et al.ll2010 
for details) and the Two Micro All Sky Survey (2MASS) catalog for the outermost region, 
in order to sample the surrouding field population. The final catalog was placed onto the 
2MASS photometric and ab solute astrometric system by following the standard procedure 
used in lFerraro et al.l (|2004| ). Photometric uncertainties are of about 0.01-0.02 and 0.03-0.05 
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for SOFI and 2MASS target, respectively. The targets have been selected along the RGB 
and in the magnitude range Kq ~ 12.3-14.0, in order to reach a suffi cient SNR (> 30- 40) 



and excluding stars brighter than the magnitude level {Kq ~12.3, see ICioni et al.l l2006l ) of 
the RGB Tip of the intermediate-age LMC population. 

Observations consists of a series of 45- minute- long exposures obtained with HRll (5600- 
5840 A and R = 24200) and HR13 (6120-6400 A and R = 22500) gratings by using the GI- 
RAFFE/MEDUSA configurations. The data reduction was performed with the GIRAFFE 
ESO pipeline that includes bias subtraction, fiat fielding, wavelength calibration and spec- 
trum extraction. Individual stellar spectra have been cleaned from the sky contribution, by 
subtracting the corresponding median sky spectrum. Finally, multiple spectra of each target 
have been coadded, reaching a SNR per pixel of ~ 40 in the faintest stars and up to ~ 100 
in the brightest ones. Identification number, right ascension and declination of each target 
are listed in Tabled! 



3. Kinematic and chemical analysis 



Radial velocities {vr) have been measured by means of the DAOSPEC code (IStetson fc Pancino 



20081 ). For each star, the spectra from the two different gratings have been analysed indepen- 
dently and the derived Vr averaged together by using the individual uncertainty as weight. 
Typical internal errors (computed as cr/ ^/Nunls 
the heliocentric correction to each Vr- 



are ~ 0.15-0.20 km s ^. Finally, we applied 



The chemic al analysis has been performed by using the suite of codes developed by 



R.L.Kurucz (see ISbordone et al.l l2004l ) aimed to computing abundances from the observed 



equivalent widths (EWs) and synthetic spectra. The model atmospheres were computed with 
the ATLAS9 code, assuming plane-parallel geometry, local thermodynamical equilibrium for 
all species and without the inclusion of the approximate overshooting in the computation of 
the convective flux. The ATLAS9 model atmosphere s were calculated with the new set of 
Opacity Distribution Functions by lCastelli fc Kuruczl (120041 ). 

The line list has been selected starting from the Kurucz/Castelli data0 and updated 
with recent laboratory data from VALD and NIST databases. We included only transitions 
with theoretical/laboratory atomic data and checked against spectral blendings through 
the inspection of suitable synthetic spectra convolved at the GIRAFFE resolution. In 
particular, we adopted for the Fe lines the atomic data from the critical compilations by 



^http:/ /wwwuser. oat. ts.astro.it/castelli/linelists. html 



- 5 - 



Fuhr. Martin &: Wiesd (Il988[ ) and iFuhr fc Wiesd (|2006[ ). that represent the most updated 
dataset of theoretical/laboratory log gf for the iron hnes. 

For some unblended transitions for which theoretical/laboratory oscillator strengths are 
not available (or for those lines not well reproduced in the solar spectrum) we derived as- 
tr ophysical oscilla t or str engths (labelled as SUN in Table [2]) by using the solar flux spectra 
of Neckel fc Labs ( 1984 ) and the model atmosphere f or the Sun computed by F. Castelljfl 
adopting the solar abundances of iGrevesse fc Sauvall (119981 ) . We estimated for these oscil- 
lator strengths an accuracy < 15%, according to the uncertainties in the line profile fitting 
and in the conti nuum placement. We decide t o inc lude in our fina l linel ist some Fe I tran- 
sitions hsted by iFuhr. Martin &: Wiesd (Il988[ ) and iFuhr fc Wiesd (120061 ) with high quoted 
u ncertainties afte r the v erification that these lines are well reproduced in the solar spectrum 
of iNeckel fc (119841) an d providing an iron abundance within ±0.1 dex from the value 
of 



Grevesse fc Sauvall (119981 ). 



Concerning th e van der Waals damping cons t ants, we adopted, whenever possible, the 
damping values by iBarklem. Piskunov fc O'Maral (120001) , while f or the other transitions the 



van der Waals constants were calculated according to iCastellil ( 120051 ). The complete line 
list used is available in Table |2l including wavelength, element code, oscillator strength (and 
corresponding accuracy), lower excitational potential and reference source. 



The atmospheric parameters have been derived as follows: 



(1) A preliminary estimate of Tes for each star was obtained from photom etric (J 
and (J — K)q colors adopti ng an average color exce ss of E{B — V) = 0.12 (jPersson et al. 



19831). the extinction la w by iRieke fc LebofskyI (119851 ) and the color-temperature calibrations 
by lAlonso et al.l (119991 ) based on the Infra- Red Flux Method. The 2MA SS magnitudes of our 
targets have been converted in the TCS photometric system (w here the lAlonso et al.l (119991 ) 
relations are defined) by using the transformations provided by ICarpentei J200lh . We used 
the photometric Teg as first-guess values and then we refined them by imposing the lack 
of any trend between Fe I abundances and the excitation potential x- Basically, we find a 
good agreement between the photometric and spectroscopic T^s, with an average difference 
of T^q'^'^-T^^°^= 139 K (a = 173 K). Finally, the spectroscopic T^g have been adopted for 
the following analysis. Typical uncertainties in the spectroscopic temperatures (calculated 
according to the slope in the plane A(Fe I) vs x) ranging from ~80 up to ~130 K. 

(2) Microturbulent velocities have been derived by requiring the lack of any trend be- 
tween Fe I abundances and the reduced EW (defined as lg{EW/X). Typical uncertainty in 



■^http: / /wwwuser. oat. ts.astro.it / castelli / sun/ap00t5777g44377klodfnew.dat 
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the Vt determination is 0.1-0.2 km/s. 



(3) Surface gravities have been estimated from the photometry, because the classical 
method to infer log g from the comparison between Fe I and Fe II lines cannot be applied, due 
to the small number of available Fe II lines in our spectra. Surface gravities have been derived 
from the Stefan-Boltzmann relation, using for each target the correspondin g spectroscopic 
Tefj. We adopted the distance modulus (m — M)q = 18.50 ( lAlved l2004j) and a typical 



evolu tionary mass (1.5 Mq) for all the targets derived from a BaSTI ( jPietrinferni et al. 
20041 ) isochrone with age of 2 Gyr, Z= 0.008 and solar-scaled chemi cal mixture, according 



to t he values es t imate d for the bulk of the LMC disk population by iBekki fc Chibal (|2005[ ) 
and ICole et al.l (|2005[ ) . The l uminosities are derive d from the Mi^-band magnitude using 



the bolometric corrections of iBuzzoni et al.l (j2010t). Internal errors in the derived log g 



(calculated by taking into account the uncertainties in the spectroscopic Teg, mass and 
luminosity) are of about 0.1-0.15 dex. The photometric estimates of temperatures and 
gravities together with those obtained from spectroscopy are listed in Table [31 

EWs were measured with DAOSPEC. The abundances of Fe and Si have been obtained 
from the EWs. O abundances have been derived from the forbidden line at 6300.31 A by 
employing spectral synthesis technique. The abundances have been determined from x^- 
minimization between the observed spectrum and a grid of synthetic spectra computed with 
the corresponding atmospheric parameters and convolved at the GIRAFFE resolution. Mg 
abundances have been derived with the same technique by using the strong line 5711.09 A; 
we decide to exclude the Mg triplet at ~ 6318 A because these weak lines are located on the 
red wing of a close auto-ionization Ca line. This produce a decrease of the continuum level 
which was not easily measurable with DAOSPEC. The final abundances are expressed in the 
usual spectroscopic fo rmalism (logarithrn i c scale in unit of solar abundances) and adopting 
the solar values from iGrevesse fc Sauvall (119981 ) except for oxygen, for which we used the 
value by ICaffau et al.l ( l2010l ) . All the abundance ratios of our sample are listed in Table ID 



4. Abundance uncertainties 

For each abundance ratio we take into account two source of uncertainties, the inter- 
nal error relative to the EW measurements and that arising from the uncertainties in the 
atmospheric parameters. 

Uncertainties for each abundance ratio due to the uncertainties in the EW measurements 
have been estimated (for each star) by dividing the line-to-line scatter by the square root 
of the number of used lines. DAOSPEC provides also for each measured spectral line an 
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error computed from the standard deviation of the local flux residuals (representing a 68% 
confidence interval on the derived value of the EW). These uncertainties have been used only 
in the computation of the slopes in the A(Fe I) vs x ^"^^ A(Fe I) vs EWR planes (see Section 
3). 



Recently, IVenn et al.l (120121 ) pointed out that the internal EW errors provided by DAOSPEC 



are under- e stima ted of a factor of 2 with respect to those derived with the classical formula 



by ICayrell (119881 ). probably due to the pixel correlation during the wavelength calibration. 
We checked this effect in our data, assuming the FWHM estimated by DAOSPEC and 
the proper SNR for each spectrum. On average, we found a general agreement between 
the two sets of uncertainties (within ± 10-20%) but not a systematic effect. However, the 
DAOSPEC errors provide an internal ranking among the measured lines in a given star and 
possible systematic under/over-estimates do not afflict our results. 

For O and Mg, for which we derived the abundances from spectral synthesis of only 
one line each, the uncertainty in the measurement has been estimated by resorting to Monte 
Carlo simulations. We study some stars as representative of the different SNR of our targets. 
We injected Poisson noise into the best-fit synthetic spectrum according to the SNR of the 
observed spectrum and we repeated the fit procedure; for each line we performed a total of 500 
Monte Carlo events, assuming as uncertainty in the fitting procedure la of the abundances 
distribution. For the O line (falling in the HR13 grating), the derived uncertainties are of 
0.04, 0.03 and 0.02 dex at SNR= 30, 50 and 80 respectively, while for the Mg line (falling 
in the HRll grating) we obtained uncertainties of 0.13, 0.09 and 0.04 dex for SNR= 30, 40 
and 70. 

We estimated also the uncertainty due to the continuum location, by using the root 
mean square of the residual provided by DAOSPEC, after removing all the fitted lines from 
the observed spectrum. For our stars the typical dispersion of the residual spectra is of 
± 3-4%, corresponding to an error in the abundances less than 0.1 dex. Keeping in mind 
that DAOSPEC estimates a global continuum along the entire spectrum (and not a local 
continuum for each individual transition), this uncertainty can be considered as a systematic 
errorbar for all the lines in a given spectrum. 

The uncertainty arisi ng from the a t mosp heric parameters has been computed following 



the approach described by lCayrel et al.l (l2004j ). The usual method to derive the errors due to 



the stellar parameters is to vary one only parameter each time, keeping the other ones fixed 
and finally adding in quadrature the derived abundance variations. Obviously, this method 
does not take into account the correlations among the atmospheric parameters, providing 
only a conservative value for the uncertainty. In our case. Teg and Vt have been optimized 
spectroscopically (and log g has been derived according to the spectroscopic temperature). 



- 8 - 



thus all the parameters are not independent each other. For each star, the temperature has 
been varied by =1= Icrr^g, because th e uncertainty in T^s dominates the derived abundances as 
pointed out by lCayrel et al.l (l2004j ). Then, we repeated the optimization procedure described 
above by keeping the temperature fixed and deriving new values for log g and Vf. The 
advantage of this m ethod is to naturally take into account the covariance terms among the 
parameters (see also IShetrone et al.ll2003l ). The procedure has been repeated independently 
for each star by considering the corresponding T^fj uncertainty computed from the error of 
the slope in the A(Fe) — x plane. Table Olists the results of this procedure for a representative 
star of our sample: the second column shows the final uncertainty for each abundance 
ratio, while the other columns are the results obtained with the usual approach to vary 
independently each parameter. Last column of Table [5] is the sum in quadrature of the 
un certainties du e to on e only parameter: these values are larger than those obtained with 
the lCayrel et al.l (120041 ) approach because the covariances among the parameters is neglected. 
Typical abundance uncertainties due to the atmospheric parameters are of the order of [Fe/H] 
= ± 0.03 dex, [0/Fe] = ± 0.04 dex, [Mg/Fe] = ± 0.06 dex and [Si/Fe] = ± 0.09 dex. 

Finally, the total internal error for each abundance ratio was obtained by adding in 
quadrature the error associated to EW measurements and atmospheric parameters. For 
[Fe/H] it turns out to be of the order of ~ 0.04-0.05 dex because of the large number of 
measured lines. 



5. Results 



The heliocentric radial velocity distribution of the stars in our sample is shown in the 
left panel of Figure [U Stars with radial velocity in the ra nge 170 < Vr < 380 km s~^ are 



considered LMC members, according to IZhao et al.l (120031 ) : only 2 stars with Vr ~100 km 
s~^ have been excluded (likely belonging to the Galaxy). The mean velocity of the sample 
is Vr= 259.3 km s^^ {a^ = 33.9 km s~^), in good a greement with previo us m easurements in 



other samples of LMC disk by I Cole et al.l (120051 ) . ICarrera et al.l ( 120081 ) and iPompeia et al. 
J2008h . 



The metallicity distribution of our sample is shown in the right panel of Figure [T] The 
entire sample has an average [Fe/H] = -0.58 dex (cxpe/H] = 0.25 dex). Also, two main 
components can be distinguished: 



1. A principal component (hereafter LMC-R) comprising ~ 84% of the observed sample 
with [Fe/H] > -0.7 dex, peaked at [Fe/H]= -0.48 dex and with a quite small dispersion 
(qpe/H] =0.13 dex). 
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2. A secondary component (hereafter LMC-P), comprising ~ 16% of the sample, peaked 
at [Fe/H] = -1.06 dex (crpc/H] = 0.18 dex) with an extended tail up to [Fe/H] ~ -1.5 
dex. 



The peak of the r netallicity d i stribu tion of LMC-R (see Fig. |2l a) is in good agreement 



with those obtained bv lCole et al. 



(120051) r|Fe/H] = -0.45 dex, cjpe/H] = 0.31 dex) (see Fig. El 



b) from the central bar and ICarrera et al.l (120081 ) ([Fe/H] = -0.47 dex (crpe/H] = 0.30 dex), 
[Fe/H] = -0.50 dex (cjpe/H] = 0.44 dex) and [Fe/H] = -0.45 dex (qpe/H] = 0.31 dex) for the 
fields located at ~ 3°, 5° and 6° at north of LMC center, respectively). A distr ibution peaked 
at a sl ightly lower value ([Fe/H] = -0.75 dex, cxpe/H] = 0.23 dex) was found by lPompeia et al. 
J2008h (see Fig. [2]). 



Figure [3] shows the behaviour of O, Mg and S function of [Fe/H]. We plotted 

for comparison the abundance ratios measured i n other environments: the Galaxy (in clud- 



ing data fo r thin disk, thick di sk and halo, by IVenn et al.l |2004J: iReddv et a. 



Letarte et al. 



2010 



(ISbordone et al. 



Mucciarelli et al 



LMC disk (IPompeia et al.l 120081). the nearby dwarf spheroi dals (jShetrone et al 



Lemasle et al. 2012; Venn et al. 



2001 



2006h. the 



2003 



2OI2I ). the Sagi ttarius dwarf galaxy 



200 



j) and t he ol d and intermediate-age LMC GCs ( Johnson et al 



2008, 2009, 2010 







2006 



The [0/Fe] ratio appears sys tematically lower than those measured in the Galaxy and 
basically consistent with those by iPompeia et al.l (120081 ) . even if they provide a few of mea- 
sures for this abundance ratio. The [Mg/Fe] ratio turns out to be subsolar in the entire 
metallicity range and ever lower than the Milk y Way stars; also, our targets show lower Mg 
abundances with respect to those measured by IPompeia et al.l (120081 ) and we attribute such 
a discrepancy to the different values of oscillator strengths and damping constants for the 
Mg line at 5711 A employing in the two analysis. Finally, the [Si/Fe] ratio seems to be barely 
consistent with the Galactic stars, despite a larger dispersion. 

Figure m shows the behaviour of the average [a/Fe] ratio, obtained by averaging the 
abundances of O, Mg and Si, as a function of [Fe/H]. We compared our targets with the 
[a/Fe] ratios of the other samples calculated by averaging the available a-element abundance 
ratios; only for the old LMC GCs we excluded O and Mg, because of the intrinsic star-to-star 
variations of these elements due to the self-enrichment process occurred in the early stages 
of these clusters (at variance to the intermediate-age LMC GCs where intrinsic spread in 



*We excluded from this comparison the LMC GCs younger than 0.5 Gyr bv lMucciarelh et al.l (|201lll2012n 
because they are associated to the last episodes of star formation, while all the targets discussed here belong 
to the RGB, thus they are older than ^1-2 Gyr. 
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O and Mg content are not detected). The overall trend of the [a/Fe] ratio with the [Fe/H] 
abundance shows a decrease at increasing metallicity. The most metal-poor stars with [Fe/H] 
< -0.7 dex show an [a/Fe] ratios larger than solar value, while stars with [Fe/H] > -0.7 dex 
show [a/Fe] ratio from solar to sub-solar values. For a given [Fe/H], the [a/Fe] ratio measured 
in the LMC appears systematically lower than those measured in the Galaxy. 

A substantial agreement of our abundances was generally found with the values mea- 
sured in dSphs, in particular with th e abundance ratios ob served in the metal-poor compo- 
nent of the Sagittarius dwarf galaxy (ISbordone et al.ll2007l ): such a similarity of a-elements, 
metallicity distribut ion and fraction of metal-poor stars has been already suggested by 
Monaco et all J2005h . 



Discussion and Conclusions 



We determined kinematic and chemical properties for 89 giant stars members of the disk 
of the LMC. This sample increases significantly the number of stars analysed so far through 
high-resolution spectroscopy (the largest sample observed by iPompeia et al.l l2008l includes 
59 stars). 

The derived metallicity distribution is dominated by LMC-R, a metal-rich and narrow 
component peaked at [Fe/H] = -0.48 dex, with a secondary component LMC-P peaked at 
[Fe/H] = ~ -1.0 dex and reaching [Fe/H] = ~ -1.5 dex (right panel of Fig. [1]). As shown 
in Figure |5l the two components show similar kinematic properties, with similar average 
(fr) (265.9 km s~^ for LMC-R and 261 km s~^ for LMC-P) and with a small increase of 
the velocity dispersion decreasing the metallicity (a^ = 25 km s~^ and = 29.7 km s~^, 
respectively) . 

The kinematic and abundance distributions can be now compared to those observed in 
the stellar population of the GC NGC 1786. All the stars of our sample belong to the RGB, 
hence their age can range from ~ 1-2 Gyr up to ~ 12-13 Gyr, thus excluding very recent (last 
500 Myr) burst of star formation. At variance to the stars belonging to a stellar cluster, 
the determination of the age for field stars through isochrone fitting is a dangerous and 
uncertain technique, because the position of a RGB star in the Color-Magnitude Diagram is 
weakly sensitive to the age and highly affected by uncertainties in the color excess, evolutive 
mass and distance. Even if precise ages for each target star cannot be determined, we 



can draw the timeline of the chemical evo 
our results with the recent simulations by 



ution in this region o f the LMC, by comparing 



Smecker-Hane et all fl2002[ ). Harris fc Zaritskvl fl2009[ ) and iRubele et al.l fl2012[ ) in different 



Bekki fc Chibal (120051) and the SFH inferred by 
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regions of the LMC; unfortunately, there are no determinations of SFH in the region of our 
targets. 

(1) NGC 1786 is an old LMC GC genera ted during the first burst of SF occurred ~ 12- 
13 Gyr ago. Following iBekki fc Chibal ( l2005l ) a loose stellar halo of old stars with a velocity 
dispersion of about ~ 40 km s~^ and a broad metallicity distribution is also expected to be 
formed during that event. 

It is interesting to check whether som e of the stars observed in our sa mple could belong to 
NGC 1786 or to the LMC halo. Note that iMucciarelh et al.l feoodboioh found for NGC 1786 
{vr) = 264.3 km s^^ {cx„ = 5.7 km s'^), [Fe/H] = -1.75 dex and [a/Fe] = 0.37 dex. Figure[5] 
shows that none of the observed stars has kinematic and chemical properties compatible 
with the population of NGC 1786. In particular, the observed sample has a systematically 
higher [Fe/H] (> -1.0 dex) and a lower [a/Fe] (< 0.20 dex), thus demonstrating that even 
LMC-P did not originated in the cluster. 

This evidence could be used to put some constraints on the presence of old metal-poor 
a-enhanced stars in the sampled LMC disk field, which turns out to be less then ~ 1%. 
Our outcome provides a substantial confirmation of the results by I Cole et al.l (120051 ) . who 
found ~ 3% of th eir sample with [ Fe/H ] < -1.5 dex (but no measurements of [a/Fe] ratio are 
available) and by lPompeia et al.l (120081 ) . who found ~ 2%, although the abundance patterns 
are incompatible with old LMC GCs. Since these samples cover different regions of the 
galaxy, the measured similar fraction of metal-poor stars in these three samples seems to 
suggest that the old, metal-poor component is distributed in a quite homogeneous way along 
the LMC disk. 

(2) After the initial burst, the LMC has undergone a long period characterized by 
a continuous SF with low efficiency. The o ccurence of this quiescent p eriod in the LMC 
field is clear ly visible in the SFH provided by ISmecker-Hane et al.l ((2002) , iHarris fc Zaritsky 
(120091 ) and iRubele et al.l (120121 ) , where a prolonged phase with very low or lacking SF is 
appreciable between ~ 4-5 Gyr an d ^ 12 Gyr ago. This quiescent period is predicted also by 
the theoretical models proposed by lBekki &: Chibal (120051 ) that derived a typical SF efficiency 
in this age range of ~ 0.1 Mq yr~^. 

During this period, the LMC evolved in isolat i on, without gravit ational interactions with 
the Galaxy and the SMC JSekki fc Chibal bood : ISesla et al.ll2007h . LMC-P likely formed 
during this period, as also suggested by the lower (by ~ 0.2-0.3 dex) [a/Fe] ratio with 
respect to the Galactic values at similar metallicity, which indicate that during this long 
period mainly SNela contribute to the gas enrichment. 



If this scenario is correct, these stars were born during the period in which no GC formed 
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(the so-called Age Gap, iRich et al.ll200ll : iBekki et al.l 120041 ). hence they are unique tracers 
of the LMC chemical evolution history between ~ 12 and 3 Gyr ago. 

We also noti ce that the LMC-P fraction (~ 16%) is larger than that (< 10%) found by 



Cole et al.l ( 20051). Thi s is pr obably due to different location of the two fields. The sample 



analysed by I Cole et al.l (120051 ) is located around the central bar (at va riance with our targets 
located in the inner disk at ~ 1.8 kpc from the LMC center). Infact. iBekki fc Chibal (l2005l ) 



suggest that the central bar formed in the last 3-4 Gyr with a marginal fraction of metal-poor 
stars in its stellar content. 

Also, it is interesting to notice that the [a/Fe] ratio of this stellar component well 
resembles the mean loc us defined by the dSph star s (grey triangles in the lower panel of Fig. |U 
but see also Fig. 11 in iTolstoy. Hill &: Tosil (120091 )). suggesting similar chemical enrichment 
histories. 

(3) LMC-R formed during the relevant burst of SF occurred in the last fe w Gyr probably 
due t o the tidal capture of the SMC by the LMC. Numerical simulations (IBekki fc Chiba 
20051 ) predict that first close encounter between the two Clouds occurred ~ 3-4 Gyr ago, when 
the LMC and SMC become a gravitational bound system. This strong interaction triggered 
th e SFR up to ^ 0.4 Mfn yr ~^. The occurrence of this SF enhancement has been confirmed 



by iHarris &: ZaritskyI (120091 ) by using Color-Magnitude Diagrams of different regions of the 
LMC. However, the binary system status of Magellanic Clouds and its timescales is still a 
matter of debate since there is no u nambiguous consensu s between dynamical simulation and 
photometric studies. For instance, iRubele et al.l ( 12012| ) analysed new photometric dataset 
from the Vista Magellanic Survey, finding an enhancement of star formation with a peak at 
~ 2 Gyr and in some regions a second peak at ~ 5 Gyr an d interpreted as th e epoch of tidal 
interaction between the LMC and the Galaxy. Conversely, iBesla et al.l (120071 ) suggested that 
the LMC-SMC system is on its first close passage about the Galaxy, being entered in the 
Milky Way virial radius only in the last 1-3 Gyr. 

LMC-R distribution is very narrow {(J[Fe/H]= 0.13 dex), suggesting that the episode of 
SF associated to the first close encounter between the two Clouds has been very efficient (in 
order to produce the majority of the LMC disk stars) and fast, because the stars have not 
had time to furtherly enrich in iron. These stars show [a/Fe] abundance ratios that are close 
to the solar value, in ni ce agreement with the values measured in the intermediate-age GCs 
( iMucciarelli et al.ll2008l ) with age of 1-3 Gyr, as shown in Figure HI 



The observed [a/Fe] ratios agree with those derived in the Sagittarius dSph giants 
(asterisks in lower panel of Fig. H]), while the nearby dSphs studied s o far do not reach such 
metallicity. The most metal-rich stars discussed by iLetarte et al.l (|2010[ ) in Fornax show 



- 13 - 



[Fe/H]~-0.6 dex, consistent with the metal-poor edge of the LMC-R population. These 
evidences (the similarity with Sgr and the difference with the nearby dSphs) seem to confirm 
that the recent SFH of the LMC has been quite complex and turbulent, and the stars of 
the LMC-R population (that dominates the metallicity distribution of our sample) are the 
products of the tidal interactions occurring among LMC, SMC and our Galaxy. 

The intrinsic dispersion of each individual abundance ratio of the LMC-R population 
suggests a high degree of homogeneity, pointing out that an efficient and homogeneous mixing 
of the SN ejecta occurs. We note that the dispersion of the [a/Fe] ratios in these stars is 
larger than the range of values covered by the intermediate-age GCs. 

Thus, the scatter in [a/Fe] abundances of the LMC-R stars reflects the larger age range 
in which these stars formed. However, we flnd a general agree between the values observed 
in the fleld and GC stars in this metallicity range, suggesting in any case a common origin 
of these stars from similar burst of SF. 

It is important to recall that there is an offset between the age of the onset of the SF 
in the LMC fleld and that of the GC formation. In fact, if the SF in the fleld started about 
4-5 Gyr ago (the precise age depends on the location in the disk and the magnitude of the 
tidal int eractions with the SM C), the GC formation restarts about 2 Gyr ago, as suggested 
also by iBekki &: Chibal ( 120051 ) . This difference is conflrmed by the mea sured ages in the 
intermediate-ages GCs (see e.g. iRich et al.ll200ll : iMucciarelli et al.ll2007al Jbl) and the lack of 
globulars in the age range between 2 and 5 Gyr. Furthermore, the lower metallicity edge of 
LMC-R stars ([Fe/H] ^ -0.7 dex), if compared to that of GCs, seems to support that the 
SF of fleld stars becomes efficient slightly before. 



Recentlv . iTsuiimoto fc Bekkil ( 12012| ) discussed the unusual low [Mg/Fe] (-0.9 dex) ratio 
measured by IColucci et al.l (120121 ) in the intermediate- age, metal-rich GC NGC 1718, sug- 
gesting that this cluster and a fraction of LMC fleld stars were formed from the metal-poor 
gas acquired by infa ll from the SMC about 1-2 Gyr ago (the so-called Magellanic Squall, 
Bekki fc Chiball2007l ). In our sample we did not detect metal-rich stars characterized by such 
low Mg or a-depletion, pointing out that the stars formed from metal-poor gas accreted by 
the SMC are (if any) a negligible (<1%) fraction of the LMC disk. 



The authors warmly thank the anonymous referee for his/her suggestions in improving 
the paper. 
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Table 1. Identification number, right ascension, declination and J, H, K magnitudes. The 

entire table is available in the online version. 



ID 


R.A. (J2000) 


Dec. (J2000) 


J 


H 


K 








(mag) 


(mag) 


(mag) 


353 


74.8429031 


-67.7529373 


14.55 


13.80 


13.61 


1415 


74.7569885 


-67.7277374 


14.52 


13.75 


13.57 


1593 


74.7403107 


-67.7479858 


14.34 


13.56 


13.36 


1954 


74.6953583 


-67.7596817 


14.11 


13.28 


13.12 


1980 


74.6899948 


-67.7497635 


14.52 


13.74 


13.58 


2353 


74.7275314 


-67.7423553 


12.56 


11.87 


11.71 


2363 


74.7172394 


-67.7453918 


14.67 


13.97 


13.82 


2376 


74.7086487 


-67.7528458 


14.62 


13.93 


13.81 


90012 


74.2391357 


-67.7102280 


13.91 


13.05 


12.93 


90119 


74.5899658 


-67.9048157 


14.27 


13.38 


13.23 



- 18 - 

Table 2. Wavelength, element, oscillator strength, excitation potential and reference 

source of adopted line list. 



A 


El. log gf accuracy E.P. 


Ref. 


(A) 


(eV) 





5586.756 


Fel 


-0.144 


<10% 


5607.664 


Fel 


-2.270 


<50% 


5611.356 


Fel 


-2.990 


<50% 


5618.632 


Fel 


-1.276 


<18% 


5619.225 


Fel 


-3.270 


<50% 


5619.595 


Fel 


-1.670 


>50% 


5624.022 


Fel 


-1.140 


<15% 


5633.946 


Fel 


-0.320 


>50% 


5636.696 


Fel 


-2.560 


>50% 


5638.262 


Fel 


-0.840 


<50% 


5640.307 


Fel 


-1.700 


<15% 


5646.684 


Fel 


-2.500 


<50% 


5650.690 


Fel 


-0.960 


>50% 


5651.469 


Fel 


-2.000 


<50% 


5652.318 


Fel 


-1.920 


>50% 


5653.867 


Fel 


-1.610 


>50% 


5661.021 


Fel 


-2.430 


<50% 


5661.345 


Fel 


-1.756 


<10% 


5662.516 


Fel 


-0.573 


<10% 


5677.684 


Fel 


-2.700 


<50% 


5678.601 


Fel 


-4.670 


<50% 


5680.240 


Fel 


-2.540 


>50% 


5691.497 


Fel 


-1.490 


>50% 


5693.640 


Fel 


-0.680 


<15% 


5704.733 


Fel 


-1.250 


<15% 


5705.465 


Fel 


-1.355 


<10% 


5714.551 


Fel 


-1.770 


<15% 


5720.886 


Fel 


-1.950 


<50% 


5731.762 


Fel 


-1.270 


<50% 


5732.296 


Fel 


-1.560 


<25% 


5741.848 


Fel 


-1.670 


<25% 


5760.345 


Fel 


-2.440 


>50% 


5767.972 


Fel 


-3.200 


<15% 


5775.081 


Fel 


-1.298 


<18% 


5776.224 


Fel 


-3.400 


<15% 


5778.453 


Fel 


-3.430 


<25% 


5793.915 


Fel 


-1.660 


>50% 



3.37 Fuhr fc Wiese (200Q) 

4.15 Fuhr. Martin fc Wiese f 19881 

3.64 Fuhr. Martin fc Wiese (1988) 

4.21 Fuhr fc Wiese (2006 1 
3.69 Fuhr. Martin fc Wiese f 19881 
4.39 Fuhr fc Wiese (2QQ6) 
4.39 SUN 

4.99 Fuhr fc Wiese Q006) 

3.64 Fuhr fc Wiese (2006) 

4.22 Fuhr fc Wiese (^006) 
4.64 SUN 

4.26 Fuhr. Martin fc Wiese (1988} 

5.09 Fuhr fc Wiese Q006) 
4.47 Fuhr. Martin fc Wiese (1988) 
4.26 Fuhr fc Wiese (2006) 
4.39 Fuhr fc Wiese f2006l 

4.58 Fuhr. Martin fc Wiese (1988) 

4.28 Fuhr fc Wiese f2006l 

4.18 Fuhr fc Wiese f2006l 

4.10 Fuhr. Martin fc Wiese f 19881 
2.42 Fuhr. Martin fc Wiese (1988} 

4.19 Fuhr fc Wiese (20061 
4.30 Fuhr fc Wiese (20061 
4.96 SUN 

5.03 SUN 

4.30 Fuhr fc Wiese f2006l 

5.09 SUN 

4.55 Fuhr. Martin fc Wiese (1988} 

4.26 Fuhr fc Wiese (2006 1 

4.99 Fuhr. Martin fc Wiese (1988} 

4.26 Fuhr fc Wiese (20061 

3.64 Fuhr fc Wiese (20061 

4.29 SUN 

4.22 Fuhr fc Wiese (20061 

3.69 SUN 

2.59 Fuhr fc Wiese (20061 
4.22 Fuhr fc Wiese (20061 
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Table 2 — Continued 



A 

(A) 


El. 


log gf 


accuracy 


E.P. 

(eV) 


Ref. 


5805.757 


Fel 


-1.590 


<25% 


5.03 


Fuhr. Martin & Wiese fl988) 


5806.725 


Fel 


-1.030 


<50% 


4.61 


Fuhr & Wiese f2006) 


5811.915 


Fel 


-2.430 


<50% 


4.14 


Fuhr. Martin & Wiese f 19881 


6120.249 


Fel 


-5.970 


<7% 


0.91 


Fuhr & Wiese r2006l 


6151.618 


Fel 


-3.299 


<10% 


2.18 


Fuhr & Wiese r2006) 


6165.360 


Fel 


-1.474 


<18% 


4.14 


Fuhr & Wiese f2006) 


6187.989 


Fel 


-1.670 


>50% 


3.94 


Fuhr & Wiese (2006) 


6200.313 


Fel 


-2.437 


<18% 


2.61 


Fuhr & Wiese f2006) 


6226.734 


Fel 


-2.220 


<50% 


3.88 


Fuhr. Martin & Wiese fl988) 


6246.319 


Fel 


-0.877 


<10% 


3.60 


Fuhr & Wiese r2006) 


6252.555 


Fel 


-1.687 


<10% 


2.40 


Fuhr & Wiese r2006) 


6322.686 


Fel 


-2.426 


<18% 


2.59 


Fuhr & Wiese f2006l 


6330.849 


Fel 


-1.190 


<15% 


4.73 


SUN 


6335.330 


Fel 


-2.177 


<18% 


2.20 


Fuhr & Wiese r2006) 


6336.824 


Fel 


-0.856 


<10% 


3.69 


Fuhr & Wiese r2006) 


6380.743 


Fel 


-1.376 


<25% 


4.19 


Fuhr & Wiese (2006) 


6300.304 


01 


-9.717 


<2% 


0.00 


Storev & ZeiDuen (2000) 


5711.095 


Mgl 


-1.724 


<10% 


4.35 


NIST 


5665.555 


Sil 


-2.040 


20% 


4.92 


Garz (1973} 


5666.677 


Sil 


-1.710 


<15% 


5.62 


SUN 


5690.425 


Sil 


-1.870 


20% 


4.93 


Garz (1973) 


5701.104 


Sil 


-2.050 


20% 


4.93 


Garz (1973) 


5793.073 


Sil 


-2.060 


20% 


4.93 


Garz (1973) 


6155.144 


Sil 


-0.870 


<15% 


5.62 


SUN 


6237.319 


Sil 


-1.100 


<15% 


5.61 


SUN 
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Table 3. Identification number, photometric temperature and gravity, radial velocity and 
spectroscopic atmospheric parameters for the stars in our sample. The entire table is 

available in the online version. 

ID TPj^°* log(7P^°' Vr logg'P''' [A/H] v^' 

(K) (dex) (kms-i) (K) (dex) (dex) (km s"!) 



353 


4074 


1.12 


269.7 


1415 


4043 


1.10 


223.6 


1593 


3989 


1.01 


254.2 


1954 


3922 


0.91 


223.3 


1980 


4018 


1.10 


272.3 


2353 


4255 


0.39 


299.4 


2363 


4227 


1.23 


265.1 


2376 


4308 


1.24 


230.8 


90012 


3906 


0.83 


243.5 


90119 


3795 


0.93 


251.1 



4150 


1.10 


-0.50 


1.50 


4200 


1.10 


-0.50 


1.80 


4150 


1.10 


-0.50 


1.60 


3850 


0.90 


0.00 


1.50 


4100 


1.10 


-0.50 


1.60 


4350 


0.40 


-0.50 


1.80 


4350 


1.30 


-0.50 


1.70 


4600 


1.30 


-1.00 


2.10 


4000 


0.80 


-1.00 


2.00 


4450 


1.20 


-0.50 


2.10 




Fig. 1. — The heliocentric radial velocity distribution (left panel) and the metallicity dis- 
tribution (right panel) of our sample of LMC dis k stars. The shaded histogra ms show the 
distribution of the 7 member stars of NGC 1786 flMucciarelh et al.lbood . boiol ). 



Table 4. Derived abundances for the stars of our sample. The entire table is available in 

the online version. 



ID 


[Fe/H] 


[0/Fe] 


[Mg/Fe] 


[Si/Fe] 




(dex) 


(dex) 


(dex) 


(dex) 



353 


-0.42 


± 





03 





01 


± 





05 


-0.16 


± 


0.08 


-0.07 


± 


0.06 


1415 


-0.62 


± 





03 





19 


± 





04 


-0.21 


± 


0.02 


0.16 


± 


0.03 


1593 


-0.58 


± 





02 





04 


± 





06 


0.06 


± 


0.09 


0.14 


± 


0.05 


1954 


-0.45 


± 





03 


-0 


18 


± 





09 


-0.13 


± 


0.10 


0.17 


± 


0.05 


1980 


-0.46 


± 





02 





06 


± 





04 


-0.22 


± 


0.10 


0.13 


± 


0.02 


2353 


-0.37 


± 





02 


-0 


38 


± 





08 


-0.16 


± 


0.09 


-0.09 


± 


0.03 


2363 


-0.50 


± 





03 


-0 


03 


± 





07 


-0.12 


± 


0.10 


0.03 


± 


0.04 


2376 


-1.08 


± 





03 





28 


± 





05 


-0.14 


± 


0.11 


0.18 


± 


0.06 


90012 


-1.24 


± 





02 





38 


± 





09 


0.02 


± 


0.07 


0.37 


± 


0.07 


90119 


-0.64 


± 





02 





24 


± 





09 


-0.11 


± 


0.13 


-0.04 


± 


0.05 



Note. — The adopted solar values are 7.50, 8.76, 7.58 and 7.55 for 
Fe, O, Mg and Si. 



Table 5. Abundance uncertainties for the star 91969. The second columri is the total 
uncertainty estimated by following the prescriptions by ICayrel et al.l (120041 ). The other 
columns show the variations in abundance due to the variation of one only parameter, 
while the last column is the sum in quadrature of these terms, without taking into account 

the covariance terms. 



Element 


Parameters 




6\og g 


6vt 


Quadrature 




uncertainty 


±100K 


±0.1 


±0.1 km s-i 






(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


Fe 


±0.01 


±0.03 


±0.02 


±0.04 


±0.05 





±0.02 


±0.01 


±0.05 


±0.02 


±0.05 


Mg 


±0.04 


±0.03 


±0.01 


±0.02 


±0.04 


Si 


±0.06 


±0.07 


±0.03 


±0.02 


±0.08 
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Fig. 2. — Comp arison o f our normaliz e d me tallicity distribution (a) with that one by 
Cole et al.l bo03( ] (b) and lPompeia et al.l fcoosh (c). 
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Fig. 3. — Behaviour of the [0/Fe], [Mg/Fe] and [Si/Fe] ratio as a functio n of [Fe/H]. Black 



dots represent our targets, the small gray diamonds are th e Galaxy data fVenn et al.l 12004 



Reddy et al.ll2006l ). the black pentagons are the LMC disk (jPompeia et al 



white squares and diamonds are t he old (IJohnson et al 



and intermediate- age LMC GCs (IMucciarelli et al. 



2008 



2006 



20081). the edged 



MucciareUi et al.l l2009l . 120101 ) 



), respectively. The [Mg/Fe] abun- 



dance ratios from iPompeia et al. f l2008h were corrected by a factor of ~ -0.11 dex to take 
into account the effect of different oscillator strength for the line at 5711 A used in both 
analysis. The average [ 0/Fe] and [Mg/Fe] ratios of the three old and metal poor GCs from 



MucciareUi et al.l (|2009l ) was obtained averaging the abundances of the stars with greater val- 



ues only, in order to avoid the effects of ant icorrelat ions. In this case the tiny black arrows 
indicate them as an "upper limit". The dark-gray arrow marks the position of NGC 1786. 
Dashed lines mark the solar value. The error bars in the bottom right corner indicate the 
cumulative (EWs + atmospheric parameters) uncertainties. 
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Fig. 4. — Behaviour of the average [a/Fe] ratio as a function of [Fe/H]. Same symbols of 
Figure [31 At variance with the intermediate-age GCs, the [a/Fe] ratio of the three old 
and metal poor GCs includes [Si/Fe] only, in order to avoid the effect due to the intrinsic 
dispersion observed in the O and Mg abundances. Dashed lines mark the solar value. The 
error bars in the top right corner indicates the average uncertainty of iron and "a-elements" 
abundances. Upper panel shows the comparison with the Galactic stars, while bottom 
panel shows the comparison with stars in the Sagittarius dwarf galaxy (aster isks) and in the 



nearby dwa rf galaxies (grey trian gles; the plotted samples include the data by lShetrone et al. 



(120011 ) an dlShetrone et al.l (120031) for Draco, Se xtant, Ursa Minor, S culpt or, Fornax, Carina 
and Leo I, Letarte et al.~(j2010 ) for Fornax and Lemasle et al. ( 2012 ) and Venn et al. (2012) 
for Carina). 



- 25 - 



-i 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 



• t 



I ■ - I 



_l I I 1 I I I I I I I 1 I I l_ 

oooooooo 

COCOCOC\JC\IC\JC\JC\J 



(^_S UJ>|) V 



Fig. 5. — The distribution of our targets (gray p oints) in the Vr - [ Fe/H| plane with of the 
7 member stars of NGC 1786 (black triangles) by iMucciarelli et al.l (120091 . 120101 ) . The black 
squares represent the average radial velocities and metallicities of the two components. 



